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The Universe is modeled as consisting of pressureless baryonic matter and a bulk vis- 
cous fluid which is supposed to represent a unified description of the dark sector. In the 
homogeneous and isotropic background the total energy density of this mixture behaves 
as a generalized Chaplygin gas. The perturbations of this energy density are intrinsically 
nonadiabatic and source relative entropy perturbations. The resulting baryonic matter 
power spectrum is shown to be compatible with the 2dFGRS and SDSS (DR7) data. 
A joint statistical analysis, using also Hubble-function and supernovae la data, shows 
that, different from other studies, there exists a maximum in the probability distribution 
for a negative present value of the deceleration parameter. Moreover, the unified model 
presented here favors a matter content that is of the order of the baryonic matter abun- 
dance suggested by big-bang nucleosynthesis. A problem of simple bulk viscous models, 
however, is the behavior of the gravitational potential and the reproduction of the CMB 
power spectrum. 
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1. Introduction 

According to the prevailing interpretation, our Universe is dynamically dominated 
by a cosmological constant A (or a dynamical equivalent, called dark energy (DE)) 
which contributes more than 70% to the total cosmic energy budget. More than 20% 
are contributed by cold dark matter (CDM) and only about 5% are in the form of 
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conventional, baryonic matter. Because of the cosmological constant problem in its 
different facets, including the coincidence problem, a great deal of work was devoted 
to alternative approaches in which a similar dynamics as that of the ACDM model is 
reproduced with a time varying cosmological term, i.e., the cosmological constant is 
replaced by a dynamical quantity. Both dark matter (DM) and DE manifest them- 
selves so far only through their gravitational interaction. This provides a motivation 
for approaches in which DM and DE appear as different manifestation of one single 
dark-sector component. The Chaplygin-gas modeffand its generalizations^ realize 
this idea. Unified models of the dark sector of this type are attractive since one 
and the same component behaves as pressureless matter at high redshifts and as a 
cosmological constant in the long time limit. While the homogeneous and isotropic 
background dynamics for the (generalized) Chaplygin gas (GCG) is well compat- 
ible with the data, the study of the perturbation dynamics resulted in problems 
which apparently ruled out all Chaplygin-gas type models except those that are ob- 
servationally almost indistinguishable from the ACDM model. ^ To circumvent this 
problem, nonadiabatic perturbations were postulated and designed in a way to make 
the effective sound speed vanish. But this amounts to an ad hoc procedure which 
leaves open the physical origin of nonadiabatic perturbations. There exists, however, 
a different type of unified models of the dark sector, namely viscous cosmological 
models, which are intrinsically nonadiabatic W In the homogeneous and isotropic 
background a one-component viscous fluid shares the same dynamics as a GCG.^121 
Now, what is observed in the redshift surveys is not the spectrum of the dark-matter 
distribution but the baryonic matter spectrum. Including a baryon component into 
the perturbation dynamics for a universe with a Chaplygin-gas dark sector, there 
appears the new problem that the unified Chaplygin-gas scenario itself is disfavored 
by the data. It is only if the unified scenario with a fixed pressureless (supposedly) 
baryonic matter fraction of about 0.043 (according to the results from WMAP and 
primordial nucleosynthesis) is imposed on the dynamics, that consistency with the 
data is obtained!^ If the pressureless matter fraction is left free, its best-fit value 
is much larger than the baryonic fraction. In fact, it becomes even close to unity, 
leaving only a small percentage for the Chaplygin gas, thus invalidating the entire 
scenario. In other words, a Chaplygin-gas-based unified model of the dark sector 
is difficult to reconcile with observations. One may ask now, whether the status of 
unified models can again be remedied by replacing the Chaplygin gas by a viscous 
fluid. It is this question that we have investigated in Refs. |4]and|8l which are sum- 
marized in this contribution. We studied the cosmological perturbation dynamics 
for a two-component model of baryons and a viscous fluid, where the latter repre- 
sents a one-component description of the dark sector. We could show that such type 
of unified model is not only consistent for a fixed fraction of the baryons but also 
for the case that the matter fraction is left free. Our analysis demonstrates that the 
statistically preferred value for the abundance of pressureless matter is compatible 
with the mentioned baryon fraction 0.043 that follows from the synthesis of light 
elements. In addition we discuss here problems of the viscous model to correctly 
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reproduce the anisotropy spectrum of the cosmic microwave backgroun and give 
a brief outlook on a possible solution in the context of causal transport theory!^ 

2. The Two-Component Model 

The cosmic medium is assumed to be describable by an energy-momentum tensor 
T tk which splits into a matter part T %k f and viscous fluid part T$ , 

T ik = pu l u k + p (g ik + , T ik = T$ + Ty, (1) 

with 

Tfl = PMU l M U k M + p M (g lk + u\ lU k M ) , Ty = P yUyUy + p V (g lk + UyUy) , (2) 

where the subscript "M" stands for matter and the subscript "V" stands for viscous. 
The total cosmic fluid is characterized by a four velocity u m while u % M represents the 
four velocity of the matter part and u v represents the four velocity of the viscous 
fluid. Energy-momentum conservation is supposed to hold separately for each of the 
components, 

1^ = 17^=0 / '/; . (3) 

Up to first order in the perturbations we have p = pu + Pv and p = Pm + Pv- 
In general, the four velocities of the components are different. We shall assume, 
however, that they coincide in the homogeneous and isotropic zeroth order, 



(background) . (4) 



Difference will be important only at the perturbative level. Let the matter be pres- 
sureless, i.e. pu = and the viscous fluid, according to Eckart's theory, be charac- 
terized by a bulk viscous pressure pv = P = — C©j where £ = const and 8 = u? i is 
the fluid expansion. Under this condition the total pressure coincides with the pres- 
sure of the viscous component. In terms of the present value qo of the deceleration 
parameter the Hubble rate can be written aJ^ 



H _ 1 



l-2 go + 2(l + g )a^ , (5) 



which coincides with the Hubble rate of a specific (a = —1/2) GCG with the general 
equation of state p = —A/p a . Since pu — Pmou 3 , we have pv = p — pmo» _3 - It is 
the total energy density that behaves as a GCG, not the component V. This type of 
unified model differs from unified models in which the total energy density is the sum 
of a GCG and a baryon component. Only if the baryon component is ignored, both 
descriptions coincide. Consequently, in the homogeneous and isotropic background, 
a generalized Chaplygin gas with a — —1/2 can be seen as a unified description of 
the cosmic medium, consisting of a separately conserved matter component and a 
bulk viscous fluid with £ = const, where the latter itself represents a unified model 
of the dark sector. 
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3. Perturbation Dynamics 

3.1. Total energy density perturbations 

The nonadiabaticity of the system as a whole is characterized by 

" P P = 3H^(--t), (6) 



p+p pp+p p \p e 

where a caret denotes a perturbation quantity. The expression ^ is governed by 
the dynamics of the total energy-density perturbation p and by the perturbations 
of the expansion scalar, which is also a quantity that characterizes the system 
as a whole. The behavior of these quantities is described by the energy-momentum 
conservation for the entire system and by the Raychaudhuri equation, respectively. 
Both of these equations are coupled to each other. The remarkable point is that 
these quantities and, consequently, the total energy density perturbation, are inde- 
pendent of the two-component structure of the medium. The reason is the direct 
relation p = — £0 between the pressure perturbations and the perturbations of the 
expansion scalar. It is convenient to describe the perturbation dynamics in terms 
of gauge invariant quantities which represent perturbations on comoving hypersur- 
faces, indicated by a superscript c. These are defined as (v is the velocity potential, 
defined by u M — v^) 

— = -+v , 5 =—. (7) 
9 9 9 

For the fractional quantities we introduce the abbreviations 

D c = , P c = . (8) 

p+p p+p 

In our case we have 

p : = t ? * ( 9) 

p e p e 

In terms of the comoving quantities the total energy and momentum balances may 
be combined into (cf. Ref. H]) 

D c -3H- D c + e c = . (10) 
9 

The expansion scalar <d is governed by the Raychaudhuri equation, which, in linear 
order, can be written in the form 



6 C + 2HO c + \AP C + — H 2 D c = . (11) 
a 2 2 

It is through the Raychaudhuri equation that the pressure gradient comes into play: 

-P — , => P c = — V -D c - — L . 
7/9 0' 27 p 2 Z^/pH 



P C = ±Z—, => P c = — ^D C -^TTeD c , (12) 
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where 7 = 1 + - . The pressure perturbation consists of a term which is proportional 
to the total energy-density perturbations D c (notice that the factor in front of D c 
is positive), but additionally of a term proportional to the time derivative D c of D c . 
The relation between pressure perturbations P c and energy perturbations D c is no 
longer simply algebraic, equivalent to a (given) sound-speed parameter as a factor 
relating the two. The relation between them becomes part of the dynamics. In a 
sense, P c is no longer a "local" function of D c but it is a function of the derivative 
D c as welP. This is equivalent to p = p(p, p). It is only for the background pressure 
that the familiar dependence p = p(p) is retained. As already mentioned, the two- 
component structure of the medium is not relevant here. 
Introducing now 

5 = 7 £> c = P - , (13) 
P 

and changing from the variable t to a, Eqs. (flQ|) and may be combined to yield 
the second-order equation 

6" + f(a)5' + g(a)5 = 0, (14) 



where 8' = and the coefficients / and g are 



/(«)= 1 



a 



3 JP 1 p 
— 6 



<„2 1 



and 



9 {a) = ~~2 



a*- 



2 p 3 7p H 2 a 2 
3 15 p 9 p 2 lp 2 



(15) 



(16) 



2 2 p 2 p 2 7 p 2 H 2 a 2 
respectively. Equation (|14[) coincides with the corresponding equation for the one- 
component case in Ref. 0] 

In Fig. Q] the density fluctuations for the viscous model are compared with those 
of the GCG model for different values of the relevant parameters. Although identical 
in the background, both models are qualitatively very different at the perturbative 
level. The density perturbations in the bulk-viscous scenario are well behaved at 
all times. The GCG model predicts (unobserved) oscillations, as was also found 
in Ref. [12l The latter behavior was the main reason for discarding these models, 
except, possibly, for very small values of a. This unwanted property does not hold 
for our viscous model. This coincides with the results of Ref. El Both models coin- 
cide for early times, confirming our previous analytical result, that non-adiabatic 
contributions are negligible in the past, but become relevant at a later period. The 
non-adiabatic contributions are essential to avoid the mentioned unrealistic features 
of GCG models. 



3.2. Relative entropy perturbations 

The relative entropy perturbations are defined by 
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Fig. 1. Absolute values (logarithmic scale) of density fluctuations as function of the scale factor 
a for v = —1 (a = 1/2) and go = —0.5 for different scales. The values of k are k = 0.5 (top left), 
k = 0.7 (top right), k = 1 (bottom left) and k = 1.5 (bottom right), all in units of hMpc' 1 . Solid 
curves represent the bulk viscous model, dashed curves the corresponding GCG model. Notice that 
both models are always different, except at very early times. 



and obey the equation 

s vm + r ( a ) s vM + s ( a )SvM = c(a)S' + d(a)S 
with the coefficients 

1 



r(a) = 



s(a) 



3 3p Pm 



2 2 p p p v +p 



3 p Pm 
a 2 p p v +p 



4 P. 



(18) 

(19) 
(20) 
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c(a) 



P 



7 Pv + P 



and 



d(a) 



P 



27a 2 pv +P 



l- P - 



P_ 

2p 



P_ 

IP 



9H 2 a 2 



IP 



9H 2 a 2 



(21) 



(22) 



The set of equations (|T4|) and (|T8|) contains the entire perturbation dynamics of the 
system. At first, the homogeneous Eq. (fl"4")l for 6 has to be solved. Subsequently, 
once 6 is known, Eq. (|18|) determines the relative entropy perturbations. 



3.3. Baryonic perturbations 

The quantity relevant for the observations is the fractional perturbation 5m = ^j- 
of the energy density of the baryons, given by 



7M = - 



PV +P 



S 



VM 



(23) 



At early times, i.e. for small scale factors a<l, the equation (|14[) has the asymp- 
totic form 



independent of and for all scales. The solutions of (|24p are 



5(a <C 1) = c\a + C2<2 



-3/2 



(24) 



(25) 



where c\ and ci are integration constants. The nonadiabatic contributions to the 
total density perturbations are negligible at high redshifts.^ For a <C 1 the coeffi- 
cients s(a), c(a) and d(a) in p^|) become negligible and r(a) — > |. Eq. (fTgf then 
reduces to 



0. 



(a <C 1) 



(26) 



It has the solution Svm = const = 0. From the definition (TT71) we find that at high 
rcdshifts 

C PM PV / -, \ / „s 

Smv = , (a < 1) . (27) 

Pm Pv 

Consequently, both the nonadiabatic contributions to the total energy-density fluc- 
tuations and the relative entropy perturbations are negligible and we have almost 
purely adiabatic perturbations 8m = 5 at a <C 1. This allows us to relate our model 
to the ACDM model at early times. We shall use the fact that the matter power 
spectrum for the ACDM model is well fitted by the BBKS transfer function!^ In- 
tegrating the ACDM model back from today to a distant past, say z = 1.000, we 
obtain the shape of the transfer function at that moment. The spectrum determined 
in this way is then used as initial condition for our viscous model. For more details 
see Ref. M 
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4. Statistical Analysis 

To estimate the free parameters of our model we perform a Bayesian analysis and 
construct the corresponding probability distribution functions. At first we consider 
large-scale-structure data from the 2dFGR£p^and SDSS DBp^ programs. The mat- 
ter power spectrum is defined by 

P k = \6 M .k\ 2 , (28) 

where 5u,k is the Fourier component of the density contrast 5m- Generally, for a 
set of free parameters {p}, the agreement between the theoretical prediction and 
observations is assessed by minimizing the quantity 

x2(p) i E M_M, (29) 

i 1 

where Nf means the number of degrees of freedom in the analysis. The quantities 
P\ h and P° bs are the theoretical and the observed values, respectively, of the power 
spectrum and Oi denotes the error for the data point i. With the help of x 2 we then 
construct the probability density function (PDF) 

P = Be-^, (30) 

where B is a normalization constant. 

To test our model against the observed power-spectra data we consider the fol- 
lowing two situations, (i) We assume the matter component to be entirely baryonic 
with a fraction CImo = 0.043 as suggested by the WMAP data. Fixing also H = 72, 
a value favored by these data as well, the only remaining free parameter is go- This 
will provide us with information about the preferred value(s) of go for the unified 
dark-sector model. Fig. [2] shows the theoretically obtained spectrum for various val- 
ues of go together with the power-spectrum data points. To better illustrate the 
relation between the predictions of the model and the observations, two different 
normalization wave numbers, fc n = 0.034/iA/pc^ 1 and k n = 0.185/iAfpc , have 
been chosen, but our statistical results do not depend on a specific normalization, 
(ii) We leave the matter fraction free, thus admitting that the matter component 
is not only made up by the baryons. This is equivalent to allow for a separate DM 
component in addition to the contribution effectively accounted for by the viscous 
fluid. This additional freedom is used to test our unified model of the dark sector 
itself. The unified model can be regarded as favored by the data if the PDF for 
the matter fraction is large around the value that characterizes the baryon fraction. 
If, on the other hand, the PDF is largest at a substantially higher value, the uni- 
fied model has to be regarded as disfavored. The results of our statistical analysis 
are visualized in Fig. [3] The center panel shows a maximum of the PDF for go at 
g ss —0.53. According to the right panel the matter-fraction probability is highest 
for Qmo ^ 0.08 and sharply decays for Omo > 0.08. 
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PS with fc v =0.034 hMpc" 1 (2dFGRS) PS with k,,=0. 1 85 hMpc" ' (2dFGRS) 




-1.6 -1.4 -1.2 -1.0 -0.8 -1.6 -1.4 -1.2 -1.0 -0.8 

Loc,„k/h(Mpc -1 ) Log ltl k/h(Mpc _l | 



Fig. 2. Power spectra (PS) normalized at k n = 0.034/iMpc~ 1 (left panels) and at 0.185/iMpc" 1 
(right panels) for different negative values of go- The top panels compare the PS with the 2dFGRS 
data, the bottom panels with the SDSS DR7 data. 




Fig. 3. Left panel: Hubble diagram, center panel: probability distribution function (PDF) for qo , 
based on a joint analysis of matter power spectrum and SN la data. The right panel shows the 
PDF for the pressureless component. 
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5. Comments 

Before summarizing our results, two comments are in order here. The first one 
concerns a more adequate, "causal" description of the viscous fluid, the second one 
points out the possibility to mimic unified models of the dark sector in terms of 
kinematic back reactions in an averaged inhomogeneous cosmology. 

5.1. Viscous fluid and gravitational potential 

While the analysis of the matter power spectrum data favors the viscous model, the 
situation changes if the spectrum of the CMB is considered. A recent investigation 
showed, that, at least for the Eckart theory applied here, the time dependence of 
the gravitational potential differs dramatically from the corresponding dependence 
of the ACDM modelfCS This confirms an earlier analysis in Ref. [§] and severely 
weakens the status of a viscous dark-sector model. Now, it is well known, that 
Eckart 's theory suffers from causality and stability problems. Therefore, one may 
hope that a more adequate description of the dark sector on the basis of causal 
thermodynamics can cure this shortcoming. The essential difference to Eckart's 
theory is, that the algebraic relation II = — £0 which we used in our analysis, is 
replaced by a differential equation 

II + tLT = -C6 , (31) 

where the relaxation time r appears as a new parameter. This parameter is related 
to the propagation speed of viscous pulses. It could be demonstrated, that with a 
very small value of this dissipative sound speed of the order < 10~ 8 , this dissipativc 
dark-sector model cou ld indeed produce a gravitational potential similar to that of 
the ACDM modelP2 

5.2. Kinematic back reaction and unified models 

There exists a line of research which tries to explain the accelerated expansion of the 
Universe as the consequence of a back reaction of a suitably averaged inhomogeneous 
cosmology. This approach, which avoids the introduction of dark energy, relies on 
the fact that the averaged Einstein tensor on the left-hand side of the field equations 
is different from the Einstein tensor of an averaged metric. This difference results in 
kinematic back-reaction terms, which modify the standard background dynamics. 
Effective fluid models of this back reaction have been constructed which include 
the Chaplygin gasffSl This seems to provide additional motivation for the further 
investigation of unified dark sector models in a different context. 

6. Conclusions 

We have modeled the cosmic substratum at the present time as a mixture of a vis- 
cous fluid and baryons. The viscous fluid is assumed to provide a unified description 
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of the cosmological dark sector. In the homogeneous and isotropic background the 
two-component system of a bulk viscous fluid and a separately conserved baryon 
component behaves as a generalized Chaplygin gas. The total energy-density per- 
turbations, however, are intrinsically nonadiabatic and coincide with those of a 
one-component viscous fluid. The fluctuations of the baryon component are ob- 
tained from a combination of the total energy density perturbations and relative 
entropy perturbations in the two-component system where the former source the 
latter. The observed matter-power spectrum is well reproduced. There do not ap- 
pear oscillations or instabilities which have plagued adiabatic Chaplygin-gas mod- 
els. The probability distribution for the deceleration parameter has a maximum 
at q w —0.53 which partially removes the degeneracy of previous studies which, 
taken at face value, were incompatible with an accelerated expansion and thus in 
obvious tension with results for the background. With the matter fraction as a free 
parameter, our analysis also revealed that the matter fraction probability is highest 
for values smaller than roughly 8%. This is a result in favor of the unified viscous 
model. Consequently, as far as the matter power spectrum is concerned, the viscous 
model remains an option for a unified description of the dark sector. On the other 
hand, the simple viscous model, based on Eckart's theory, has problems to account 
for a gravitational potential that reproduces the CMB power spectrum. Possibly, 
this problem may be solved with the help of a causal transport theory for the bulk 
viscous pressure. 
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